Test of the Additivity Principle for Current Fluctuations 
in a Model of Heat Conduction 
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The additivity principle allows to compute the current distribution in many one-dimensional (ID) 
nonequilibrium systems. Using simulations, we confirm this conjecture in the ID Kipnis-Marchioro- 
Presutti model of heat conduction for a wide current interval. The current distribution shows both 
Gaussian and non-Gaussian regimes, and obeys the Gallavotti- Cohen fluctuation theorem. We verify 
the existence of a well-defined temperature profile associated to a given current fluctuation. This 
profile is independent of the sign of the current, and this symmetry extends to higher-order profiles 
and spatial correlations. We also show that finite-time joint fluctuations of the current and the 
profile are described by the additivity functional. These results suggest the additivity hypothesis as 
a general and powerful tool to compute current distributions in many nonequilibrium systems. 
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Nonequilibrium systems typically exhibit currents of 
different observables (e.g., mass or energy) which char- 
acterize their macroscopic behavior. Understanding 
how microscopic dynamics determine the long-time av- 
erages of these currents and their fluctuations is one of 
the main objectives of nonequilibrium statistical physics 
[H, [2, H, i, B @, Q- This problem has proven to be 
a challenging task, and up to now only few exactly- 
solvable cases are understood [H, [3, An important 
step in this direction has been the development of the 
Gallavotti-Cohen fluctuation theorem Q, which re- 
lates the probability of forward and backward currents 
reflecting the time-reversal symmetry of microscopic dy- 
namics. However, we still lack a general approach based 
on few simple principles. Recently, Bertini and coworkers 
[l| have introduced a Hydrodynamic Fluctuation Theory 
(HFT) to study large dynamic fluctuations in nonequi- 
librium steady states. This is a very general approach 
which leads to a hard optimization problem whose solu- 
tion remains challengin^in most cases. Simultaneously, 
Bodineau and Derrida [2] have conjectured an additiv- 
ity principle for current fluctuations in ID which can be 
readily applied to obtain quantitative predictions and, 
together with HFT, seems to open the door to a general 
theory for nonequilibrium systems. 

The additivity principle (also referred here as BD the- 
ory) enables one to calculate the fluctuations of the cur- 
rent in ID diffusive systems in contact with two bound- 
ary thermal baths at different temperatures, ^ Tj^. 
It is a very general conjecture of broad applicability, 
expected to hold for ID systems of classical interact- 
ing particles, both deterministic or stochastic, indepen- 
dently of the details of the interactions between the par- 
ticles or the coupling to the thermal reservoirs. The 
only requirement is that the system at hand must be 
diffusive, i.e. Fourier's law must hold. If this is the 
case, the additivity principle predicts the full current 
distribution in terms of its first two cumulants. Let 



V N{q,TL,T}i,t) be the probability of observing a time- 
integrated current Qt — qt during a long time i in a 
system of size N . This probability obeys a large devia- 
tion principle i, PAr(g, Tl, Tr, t) ~ exp[iJPjv(g, Tl, Tr)], 
where !Fn{<1,Tl^Tr) is the current large-deviation func- 
tion (LDF), meaning that current fluctuations away 
from the average are exponentially unlikely in time. 
The additivity principle relates this probability with the 
probabilities of sustaining the same current in subsys- 
tems of lengths N — n and n, and can be written as 
^N{q, TljTr) = maxT [TN-n{q, Tl,T) + Tniq, |^^-r)] 
for the LDF [2]. In the continuum limit one gets [3| 



G(q) = — min 



r-l [q + K[T,{x)]mx)y 



2a[T,(a 



■da; 



(1) 



with Q{q) — NTNij^) ~ drop the dependence on 
the baths for convenience -, a; G [0, 1], and where k{T) 
is the thermal conductivity appearing in Fourier's law, 
{Qt)/t = —k{T)VT, and a-{T) measures current fluctu- 
ations in equilibrium (Tl = Tr), {Qt)/t = <j{T)/N. The 
optimal profile Tq{x) derived from ([T]) obeys 



da; 



q^{l + 2Ka[Tq{x)]} , (2) 



where K{q^) is a constant which fixes the correct bound- 
ary conditions, T,(0) = Tl and r,(l) = Tr. Eqs. ^ and 
^ completely determine the current distribution, which 
is in general non-Gaussian and obeys the Gallavotti- 
Cohen symmetry, i.e. Q{~q) = Q{q) ^ E q with E some 
constant defined by k{T) and cr(T) 3- 

The additivity principle is better understood within 
the context of HFT [Jl j which provides a variational prin- 
ciple for the most probable (possibly time-dependent) 
profile responsible of a given current fluctuation, leading 
usually to unmanageable equations. The additivity prin- 
ciple, which on the other hand yields explicit predictions. 
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FIG. 1: (Color online) Main: ^i(A) for the ID KMP model. 
Fluctuations are Gaussian for A « 0, but non-Gaussian in the 
tails. Inset: Test of the Gallavotti-Cohen fluctuation relation. 



is equivalent within HFT to the hypothesis that the opti- 
mal profile is time- independent, an approximation which 
in some special cases breaks down for for extreme cur- 
rent fluctuations Q. Even so, the additivity principle 
correctly predicts the current LDF in a very large cur- 
rent interval, making it very appealing. 

The goal of this paper is to test the additivity prin- 
ciple in a particular system: the ID Kipnis-Marchioro- 
Presutti (KMP) model of heat conduction [9(. The model 
is defined on a ID open lattice with N sites. A config- 
uration is given by {ej, i = 1 . . . N}, where G M+ is 
the energy of site i. Dynamics is stochastic, proceeding 
through random energy exchanges between randomly- 
chosen nearest neighbors. In addition, boundary sites 
(i = 1, TV) may also exchange energy with boundary heat 
baths whose energy is randomly drawn at each step from 
a Gibbs distribution at the corresponding temperature 
(Tl for i = 1 and Tr for i = N). For Tl Tr KMP 
proved that the system has a steady state characterized 
by an average current {q) and a linear energy profile 
T.t[x) ^ Tl + x{Tr - Tl) in the iV ^ oo hydrody- 
namic scaling limit, such that Fourier's law holds with 
k(T) = i. Moreover, one also gets a{T) = T^. This 
model plays a fundamental role in nonequilibrium sta- 
tistical physics as a benchmark to test new theoretical 
advances, and represents a large class of quasi-lD dif- 
fusive systems of technological and theoretical interest. 
Furthermore, the KMP model is an optimal candidate to 
test the additivity principle because: (i) One can obtain 
explicit predictions for its current LDF, and (ii) its sim- 
ple dynamical rules allow a detailed numerical study of 
current fluctuations. 

Analytical expressions for Q{q) and Tq{x) in the KMP 
model, though lengthy, are easily derived from the ad- 
ditivity principle, see eqs. ([T]) and ([2]), and will be 
published elsewhere We find that PA,(q, T^, Tr, t) 

is Gaussian around {q) with variance cr(T), while non- 



Gaussian tails develop far from {q). Exploring by 
standard simulations these tails to check BD theory 
is unapproachable, since LDFs involve by definition 
exponentially- unlikely rare events. Recently Giardina, 
Kurchan and Peliti ll| have introduced an efficient 



method to measure LDFs in many particle systems, 
based on a modification of the dynamics so that the rare 
events responsible of the large deviation are no longer 
rare [l^. This method yields the Legendre transform 
of the current LDF, /x(A) = N'^ ma3iq[g{q) + \q]. The 
function /z(A) can be viewed as the conjugate potential to 
Q{q), with A the parameter conjugate to the current g, a 
relation equivalent to the free energy being the Legendre 
transform of the internal energy in thermodynamics, with 
the temperature as conjugate parameter to the entropy. 

We applied the method of Giardina et al to measure 
^(A) for the ID KMP model with AT = 50, = 2 and 
Tr = 1, see Fig. [T] The agreement with BD theory 
is excellent for a wide A-interval, say —0.8 < A < 0.45, 
which corresponds to a very large range of current fluc- 
tuations. Moreover, the deviations observed for extreme 
current fluctuations are due to known limitations of the 
algorithm [l^, so no violations of additivity are ob- 
served. In fact, we can use the Gallavotti-Cohen sym- 
metry, ^(A) = - E) with E = (T-i - T-i), to 
bound the range of validity of the algorithm. The inset 
to Fig. [l]shows that this symmetry holds in the large cur- 
rent interval for which the additivity principle predictions 
agree with measurements, thus confirming its validity in 
this range. However, we cannot discard the possibility of 
an additivity breakdown for extreme current fluctuations 
due to the onset of time-dependent profiles [1| , although 
we stress that such scenario is not observed here. 

The additivity principle leads to the minimization of a 
functional of the temperature profile, Tq[x), see eqs. ([T]) 
and ([2]). A relevant question is whether this optimal pro- 
file is actually observable. We naturally define Tq{x) as 
the average energy profile adopted by the system during 
a large deviation event of (long) duration t and time- 
integrated current g, measured at an intermediate time 
1 <^T -^t. Top panel in Fig. [5] shows Tq{x) measured in 
standard simulations for small current fluctuations, and 
the agreement with BD predictions is again very good. 
This confirms the idea that the system modifies its tem- 
perature profile to facilitate the deviation of the current. 

To obtain optimal profiles for larger current fluctua- 
tions we may use the method of Giardina et al [ll| . This 
method naturally yields T^'^{x), the average energy pro- 
file at the end of the large deviation event {t — t)^ which 
can be connected to the correct observable Tx{x) by 
noticing that KMP dynamics obeys the local detailed bal- 
ance condition Q , which guarantees the time reversibility 
of microscopic dynamics. This condition implies a sym- 
metry between the forward dynamics for a current fluc- 
tuation and the time-reversed dynamics for the negative 
fluctuation |T3| that can be used to derive the following 
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FIG. 2: (Color online) Excess profiles for different currents 
measured in standard simulations (top) and with the method 
of [ll[ (bottom). Much larger current fluctuations can be 
explored in the latter case. The inset shows a scaling plot of 
excess profiles for small current deviations, q G [0.44, 0.62]. In 
all cases lines correspond to BD theory. 



relation between midtime and endtime statistics fld^ 

Px{C) = Ax {p'c^)-'Pr''{C)Plf-EiC) . (3) 

Here Px'^'^iC) [resp. Pa(C)] is the probability of con- 
figuration C at the end (resp. at intermediate times) 
of a large deviation event with current-conjugate pa- 
rameter A, and p"^ ~ '^iLi Pi^i] is an effective 
weight for configuration C = {ei,i = 1...7V}, with 
Pi — + Ej^^, while A is a normalization constant. 
Eq. ([3]) imphes that configurations with a significant con- 
tribution to the average profile at intermediate times are 
those with an important probabiUstic weight at the end 
of both the large deviation event and its time-reversed 
process. Using the relation ^ and a local equilibrium 
(LE) approximation [l^l, we obtained T\(x) for a large 
interval of current fluctuations using the method of Gi- 

[2 (bottom). 



ardina et al [ll|, see Fig. [2] (bottom). The excellent 
agreement with BD profiles confirms the additivity prin- 
ciple as a powerful conjecture to compute both the cur- 
rent LDF and the associated optimal profiles. Moreover, 
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FIG. 3: (Color online) Main: Fluctuations of total energy vs 
q measured in standard simulations for t = 10® (Q) and LE 
results (□). Inset: Similar results for t — 4000. Notice the 
non-trivial structure. 



this agreement shows that corrections to LE are weak in 
the KMP model, though we show below that these small 
corrections are present and can be measured. 

An important consequence of cq. ([3]) is that P\{C) = 
P-x-e{C), or equivalently Pq{C) — P-q{C), so midtime 
statistics does not depend on the sign of the current. 
This implies in particular that Tq{x) — T^q{x), but also 
that all higher-order profiles and spatial correlations are 
independent of the current sign. 

As another test, BD theory predicts for q ~ {q) = ^ 
the limiting behavior 



Tq{x) - Tstjx) 

2q-l 



x{l - x){5 - x) + Oi2q - 1) . 



The inset to Fig. [2] confirms this scaling for Tq (x) and 
many different values of q around its average. 

We can now go beyond the additivity principle by 
studying fluctuations of the system total energy, for 
which current theoretical approaches cannot offer any 
prediction. An exact result by Bertini, Gabrielli and 
Lebowitz fl!] predicts that m2(e) = r7i2^(e) -I- ^^^(Tl - 
Tr)^ , where m2(e) is the variance of the total energy, 
7712 ^ is the variance assuming a local equilibrium (LE) 
product measure, and the last term reflects the correction 
due to the long-range correlations in the nonequilibrium 
stationary state [14|. In our case, 7773^ = (T| + T^Tr + 
r|)/3 = 7/3 « 2.3333, while 7772 = 29/12 « 2.4166. Fig. 
Oplots 77i2(e, q) = N[{e^){q) — (e)^(g)] measured in stan- 
dard simulations, showing a non-trivial, interesting struc- 
ture for 7772 (e, q) which both BD theory and HFT cannot 
explain. One might obtain a theoretical prediction for 
7712(6,(7) by supplementing the additivity principle with 
a LE hypothesis, i.e. Pq{C) « B^Ii exp[-e,/Tq( j^)], 

which results in m,2^{e,q) — J^dxTq{x)'^. However, 
Fig. [3] shows that m2^{e, (q)) « 2.33 as corresponds 
to a LE picture, and in contrast to the exact result 
7772(e, (q)) w 2.4166. This proves that, even though LE is 
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FIG. 4: (Color online) Main; Finite-time profile fluctuations 
for N = 50 and extended BD prediction. 



a sound numerical hypothesis to obtain Tx{x) from end- 
time statistics using the method of Giardina et al 13, U | , 
see Fig. [2] (bottom), corrections to LE become apparent 
at the fluctuating level. 

For long but flnite times, the profile associated to a 
given current fluctuation is subject to fluctuations itself. 
These joint fluctuations of the current and the profile are 
again not described by the additivity principle, but we 
may study them by extending the additivity conjecture. 
In this way, we now assume that the probability to flnd 
a time- integrated current q/N and a temperature profile 
Tq[x) after averaging for a long but finite time t can be 
written as pNij^ ,fq{x);t] ~ e^Y>{jqQ['i,Tq{x)]}, where 
Q[q,Tq{x)] is the functional of eq. ([T]) but not subject to 
the minimization procedure (with the notation change 



Tq ^ Tq). In this scheme the profile obeying eq. ([2]), 
i.e. the one which minimizes the functional Q, is the 
classical profile Tq{x). We can make a perturbation of 
Tq{x) around its classical value, Tq{x) — Tq{x) -\- riq(x), 
and for long t the joint probability of q and riq{x) obeys 



PN{q,t) 



exp 



1 



dx(lyAq{x,y)'qq{x)-qq{y) 



where the xy-symmetric kernel is Aq{x,y) = 
t(47Vr2)-i [2Tq-%d. -dl- %q^K{q^)\ Six - y). 
In order to check this approach, we studied the ob- 
servable {f^{x)) -T^{x) = (?7q(x)) = A-^{x,x). The 
function A~^{x,x) can be written in terms of the 
eigenvectors and eigenvalues of kernel Aq{x,y) [T^]. For 
the particular case g = (q) = ^ we were able to solve 
analytically this problem in terms of Bessel functions, 
leading to an accurate numerical evaluation of Aq^(x, x) 
[lo| . We compare these results in Fig. |4]with standard 
simulations for both i = 4 x 10^ and t — 10^. Good 
agreement is found in both cases, pointing out that BD 
functional Q[q^Tq(xy\ contains the essential information 
on the joint fluctuations of the current and profile. 



In summary, we have confirmed the additivity princi- 
ple in the ID KMP model of heat conduction for a large 
current interval, extending its validity to joint current- 
profile fiuctuations. These results strongly support the 
additivity principle as a general and powerful tool to com- 
pute current distributions in many ID nonequilibrium 
systems, opening the door to a general approach based 
on few simple principles. Our confirmation does not dis- 
card however the possible breakdown of additivity for 
extreme current fluctuations due to the onset of time- 
dependent profiles, although we stress that this scenario 
is not observed here and would affect only the far tails 
of the current distribution. In this respect it would be 
interesting to study the KMP model on a ring, for which 
a dynamic phase transition to time-dependent proflles 
is known to exist [ij. Also interesting is the possible 
extension of the additivity principle to low-dimensional 
systems with anomalous, non-diffusive transport proper- 
ties 0, or to systems with several conserved fields or in 
higher dimensions. 
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